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ABSTRACT 



We present the detection of bright [CII] emission in the z=4.76 submillimetre galaxy LESS J033229.4-275619 using the Atacama 
Pathfinder Experiment. This represents the highest redshift [CII] detection in a submm selected, star-formation dominated system. The 
AGN contributions to the [CII] and far-infrared (FIR) luminosities are small. We find an atomic mass derived from [CII] comparable 
to the molecular mass derived from CO. The ratio of the [CII], CO and FIR luminosities imply a radiation field strength Go~10^^ and 
a density ~10'* cm"'' in a kpc-scale starburst, as seen in local and high redshift starbursts. The high L[cii]/Lfir=2.4x10"' and the very 
high Z.[cii]/^'CO{1-0)~10* are reminiscent of low metallicity dwarf galaxies, suggesting that the highest redshift star-forming galaxies 
may also be characterised by lower metallicities. We discuss the implications of a reduced metallicity on studies of the gas reservoirs, 
and conclude that especially at very high redshift, [CII] may be a more powerful and reliable tracer of the interstellar matter than CO. 

Key words, galaxies: high-redshift - galaxies: ISM ~ submillimetre ~ infrared: galaxies 



1. Introduction 

Blank field submillimetre (submm) surveys have uncovered a 
population of distant galaxies undergoing a major epoch of 
star formation. To dat e, hundreds o f thes e submm g alaxies 
(SMGs) are known (e. s. 'Coppin et al.', "20061 IWeiB et all |2009|). 
However, we only have basic knowledge of the star forma- 
tion process in SMGs, and how it compares with star forma- 
tion regions in the local universe. The most obvious tracer of 
the cold molecular gas from which stars f orm is CO, which 
has been extensivel y observed in SMGs (e.g. lGreve et al.ll2005t 
iTacconi et al.L 12006) . and ha s now also been detected in other 
high redshift galaxies (e.g. iDaddi et all l2009t ITacconi et al.L 
I2OIOI) . These CO observations have shown that SMGs are 
gas rich systems {M{H2) - 10'" - IO^'Mq), providing suf- 
ficient fuel to sustain massive star formation rates exceeding 
1000 Moyr '. While the CO lines offers numerous advantages, 
notably their brightness allowing t heir emission to be spatially 
resolved (e.g. ITacconi et al.L 1200 8). there are also some com- 
plications. First, there are uncertainties of a factor five in the 
CO to H2 conversion factor Xco, especially when integrated 
over an e ntire galaxy where significant variations in Xqo are 
expected jShettv et al.L 1201 lb . Second, the derived gas masses 
depend on the excitation conditions, and a well sampled CO 
spectral line energy distribution is required to determine the 
full distribution of temperature and de nsity components (e.g. 
IWeifi et al.Ll2007UDanie lson et al.','201 iV Third, significant dust 
optical depth effects (|Papadopoulos et al., 2010) may suppress 
the high excitation lines, while contributions from shock excita- 
tion ( Papadopoulos et al.L 12008) may boost them. 

The [CII] /I157.74/zm line has recently emerged as a pow- 
erful alternative line to study the interstellar medium in high 
redshift objects. It arises from photodissociation regions (PDRs) 



associated with star-forming regions and is the dominant cool- 
ing line, representing up to 1% of the total luminosity (e.g. 
ICrawford et all \l9S^ IStacev et al.l [MB- The [CII] Hne is 
nearly always optically thin, though it may also be attenu- 
ated by dust extinction. [CI I] detections have now been re- 



porte d in five z > 2 o bjects (Maiolino et al., 2005; lono et al 
20061: iMaiohno et all i2009; Wagg et al., 2010; Iv ison et al 



20101). and 12 galax ies at 1 <z<2 dHailev-Dunsheath et alil201C ; 
Stacev et all l2010l) . The first high-redshift [CII] observations 



seemed to confirm the trend seen in local galaxies, where 
the ratio of [CII] and far-infrared (FIR) luminosity L[cii]/iFiR 
is lower by ab out an order of magnitude for sources with 
Lfir > 10"Lo (ILuhman et all Il998l llool . However, subse- 
quent observations revealed that about half of high redshift [CII] 
lines have equal relative brightness than nearby normal galax- 
ies. Recentlv. iGracia-Carpio et al.l (|201 ll) showed that z>l and 
local galaxies overlap when using Lfir/^h, ratio instead of 
Lfir as a proxy for L[cii]/Lfir, probably because Lfir/^Hj is 
more closely related to the properties of the clouds th an Lfir 
alone. The improved [CII] statistics also allowed Stace v et"an 
(120 lot) to conclude that starburst dominated systems have simi- 
lar L[cii]/Lfir than local normal galaxies, while AGN dominated 
systems have lower ratios like those in local ultra-luminous in- 
frared galaxies (ULIRGs). Both classes imply kpc-size emit- 
ting regions, but the AGN-dominated sources app ear to have 8 
times more intense far-UV radiation fields, which fStacev et all 
(2010) speculate result s from an AGN-t riggered population of 
very young stars. Mai olino et all (l2009l) suggested that lower 
metallicities may explain the increased relative [CII] brightness 
in the highest redshift sources compared to similarly bright FIR 
sources in the local universe. 

Clearly, the presence of an AGN complicates the determi- 
nation of the physical conditions in the star-forming regions 
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of their host galaxies, as the AGN may contribute to both the 
AGN s may also enrich the ISM with 



Redshift 



FIR continuum a nd [CII 

metals (e.g. Juarez et al.L[2009h . which would complicate study- 
ing th e impact of metallicity in SMGs. In a recent Herschel 
study. ISantini et all (1201 Q) found that l<z<2 SMGs have higher 
dust masses than local ULIRGs, while their derived gas metal- 
licities are lower. It is as yet unclear if this apparent dis- 
crepancy continues to the highest redshift SMGs, who are ex- 
pected to have lower metallicities. In this letter, we present 
[CII] observations of LESS J033229.4-275619atz=4.755 (here- 
inafter L ESS J033229.4; poppin et al., 2009). After AzTEC-3 
(z=5.298: 'Rieche rs et all Hoift iCapak et al.l. |20TT). this is the 
most distant SMG currently known. We detect bright [CII] emis- 
sion and show it to be consistent with an extended starburst re- 
gion, and having close characteristics to local sub-solar metal- 
licity regions. Throughout, we assume Ho = 71 km s"' Mpc"', 
Qa = 0.73, and Q,n = 0.27. We define the far-IR luminosity as 
integrated over 42. 5-122.5 //m, to al low an easy compar ison with 
iHailev-Dunsheath et al.1 (l2010l) and lStacev etaP (l2010h . 



2. Observations and Results 

LESS J033229.4 was observed with the Swedish Heterodyne 
Facility Instrument (SHFI; Vassiley et_al-L'2008) on the Atacama 
Pathfinder Experiment (APEXQ; Gusten et al., 2006) 12 m tele- 
scope in six observing runs, from 2010 April 16 to August 
24 (run ID 085.A-0584(A)), for a total of 14.5 h of observa- 
tions, including sky observations, calibrations and overheads. 
The weather conditions were generally very good with precip- 
itable water vapour 0.4<PWV<1.0mm. We used the APEX-2 
receiver tuned to 330.242937 GHz (in lower sideband), which is 
the expected frequency of the [CII]/il57.741 yum li ne at the red- 
shift z=4.762 provided by the CO(2-l) detection dCoppin et al.L 
12010.) . The two units of the Fast Fourier Transfer Spectrometer 
(EFTS) backend were offset by +200 km s ' to increase the ve- 
locity coverage. However, as the line was well centred at the 
expected position, we discarded the outer 200 km s ' covered 
only by a single EFTS as we found them to be noisier and domi- 
nated by baseline subtraction effects. Observations were done in 
wobbler-switching mode, with a symmetrical azimuthal throw 
of 20" and a frequency of 0.5 Hz. Pointing was checked on the 
nearby source R-For every ~1 h and found to be better than 3" 
(with a beam size of 20"). The focus was checked on the avail- 
able planets, especially after sunrise when the telescope defor- 
mations are largest. 

The data were analyzed by using CLASS (within the 
GI LDAS-IRAM p a ckage ), following the procedures described 
bv iMaiolino et al J (l2009h . In short, we determined the noise 
in the individual spectra using a zero order baseline, and re- 
moved the noisiest 20%. We subtracted a first order base- 
line to the combined spectrum by interpolating the channels 
in the line-free velocity ranges -850<AV'<-600kms ' and 
-250<Ay<-50kms The final spectrum includes a total of 
2.5 h on-source integration time. 

Figure [1] shows the resulting spectrum. The [CII] line is de- 
tected at the 7.6cr level. The significance of our detection is fur- 
ther strengthened by the correspondence in line width with the 
CO(2-l) detection. Table [1] compares the redshifts of the Lya, 
[Oil], CO(2-l) and [CII] lines, derived from single Gaussian 
fits. The [CII] is offset by -88±74kms-i, +130±280kms-i and 
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Fig. 1. Top panel: AT CA CO(2-l) spectru m of LESS J033229.4 
binned to 90km s^' (ICoppin et all 1201 Ol) . Bottom panel: [CII] 
spectrum obtained with SHFI/ APEX binned to 28kms ' 
(31 MHz). The blue solid lines shows the best fitting Gaussian 
profiles. The bo ttom axis shows the velocity offset relative to the 
Lyo' redshift of ICoppin et al.l (l2009l) . which is also marked as a 
dotted vertical line. 



Table 1. Obser-ved properties of LESS J033229.4. 



Parameter 


Value 


Reference 




4.762±0.002 


CoDpin et al. (2009) 


Z[OII] 


4.751±0.005 


Alaghband-Zadeh et al., in prep. 


ZCO 


4.755±0.001 


CoDDin et al. (2010) 


Z[CII] 


4.7534±0.0009 


this paper 


AVco 


160±65kms-' 


CoDDin et al. (2010) 


AV[cii] 


161±45kms-' 


this paper 


'co 


0.09±0.02 Jykms-' 


CoDDin et al. (2010) 


/[CII] 


14.7±2.2Jykms-' 


this paper 


^[CIIJ 


1.02 ±0.15 X 10'"Lo 


this paper 


^^€0(2-1) 


7.6± 1.7x IO'^Lg 


CoDpin et al. (2010) 


ft 


4.2+ 1.8x lO'-Lo 


Coppin et al. (2009^ 



^Converted using L. 



'41.5-l22.Sfim 



^S-lOOOiim/ 



' APEX is a collaboration between the Max-Planck-Institut fur 
Radioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory. 



-470± 120 km s ' from CO(2-l), [Oil] and Lya, respectively. 
The red velocity offset of t he Lyg hne ha s been observed in many 
high redshift objects (e.g. ISteidel et al.L l2010l) . and can be ex- 
plained by HI absorption in the interstellar medium surrounding 
the galaxy. The other 3 lines are more likely to trace the systemic 
redshift. Because (i) the velocity offset between [CII] and CO is 
only 1.2cr, and (ii) none of the z>l [CII] detections show signifi- 
cant offsets, we ascribe the offset in LESS J033229.4 to the poor 
S/N of the CO data rather than to a real physical velocity offset 
between the [CII] and CO emitting regions. 

3. Discussion 

The [CII] emission from LESS J033229.4 is among the brightest 
observed at z>2, in particular when compared to the CO emis- 
sion. The observed L[cn]/LFiR=2.4xlO"-' is similar to local nor- 
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mal galaxies rather than local ULIRGs. Even when considering 
iFiR/MH,=260, LESS J033229.4 still falls well above the re- 
lation defined by lo cal and high redshift objects (Figure 1 of 
iGracia-Carpioetail l20TTh . We now use the CO, [CII] and FIR 
luminosities to examine the physical conditions of the ISM in 
LESS J033229.4 and try to explain the blight [CII] emission rel- 
ative to other populations at both low and high redshift. 

The L[cii]/Lfir versus Lco(i-0)/i'FiR diagram is a particu- 
larly powerful tool as the two ratios are sensitive to gas density 
n a nd the inci dent far-ultraviolet (FUV; 6eV</iv<13. 6 eV) flux 
Go dStacev et al., 1991; Hailev- Dunsheath et aUl2010l) . Figure|2] 
compares LESS J033229.4 with other objects and the solar 
metallicity PDR model curves of lKaufmaneTaP il99% . This 
diagram can be used to roughly distinguish both the n and Go 
appropriate for a source, but with some important caveats. First, 
the CO luminosity is expressed in terms of the ground-state tran- 
sition of CO, while for LESS J033229.4 and BRI 1335-041:0, 
only the (2- 1 ) trans ition has been observed. For consistency with 
IStacev etalJ (l2010h . we assume Lco(2-i)/ico(i-0)=7.2 (90% of 
the thermalised, optically thick value) to calculate the Lco(i-O)- 
Second, Figure |2] assumes that the contributions to the [CII] 
flux from the diffuse ionised medium and from cosmic rays 
(CR) heated gas are small compared to th e ones from PDRs. 
Using PDR models, iMeiierink et al.1 (1201 ih predict that the C° 
and C^ abundances increase rapidly with the CR rates, es- 
pecially when compared to CO . If the CR rate is related to 
the star formation rate (Papado poulosI 1201 Ol) . this should thus 
be reflected in a higher [CII]/CO ratio, while powerful star- 
bursts have similar ratios than normal galaxies (see Fig. |2). 
Any contribution from CR ionisation therefore does not ap- 
pear to dominate the [CII]/CO ratio. Third, we assume that 
the AGN does not contribute significantly to the FIR and [CII] 
luminosity. LESS J033229.4 does contain a moderate lumi- 
nosity AGN, revealed by (i) ultr a-deep X-ray spe ctroscopy 
revealing a Compton-thick AGN (iGiUi et all l201lb . and (ii) 
the detection of [NV] A\240A emission dCoppin et all l2009h . 
From their SED fitting, iGilli et all (1201 ll) derive an AGN bolo- 
metric luminosity of 7x10"*^ ergs"', which represents between 
20 and 30% of the total luminosity derived between 8 and 
lOOOyum (and probably lower if the real LpiR has been under- 
estimated). It is therefore safe to assume that the FIR emis- 
sion in LESS J033229.4 is starburst dominated. The [CII] emis- 
sion could also contain a contribution from the AGN through 
the X-ray Dominated Regions (XDRs) in the molecular torus 
surrounding the AGN. IStacev et aP (1201 Oh estimate the XDR 
contributions using a scaling relation L[cn],AGN=2xlO"-'Lx-ray 
Using the intrisic (de-absorbed) re st-frame 2-lOkeV luminosity 
2.5x10'*'' erg s ' (iGilh et all l201lh . we find that the predicted 
XDR contribution is ~1.3% of the observed [CII] lumonisity. 

In our PDR model analysifl we multiply the CO flux by 
a factor of two to correct for geometry effects, and scale the 
[CII] flux by a factor 0.7 to remove non-PDR contr i bution s 
(shown by a vector in Fig.|2l iHailev-Dunsheath et al.L l2010t) . 
The geometry correction applies to all external galaxies where 
the FUV illumination comes from all sides, so this will not 
impact the relative position of LESS J033229.4 compared to 
other galaxies. We derive 10 cm and Go~10\ Such a FUV 
field is on the high end of those seen in low redshift nor- 



2 We use the M(H2)= 1.6x10'" Mq from lCoppin et al.l ( l2010h . derived 
using A'co=0.8Mo (Kkms"' pc-)-'. 

^ Th e [CII] detection of this source was published after IStacev et alJ 
( 120101) : we add the data of this target to Figure |2] 

"* We used the PDR toolbox code at http ://du stem. astro.umd. eduj] 




10"' 10"'^ 10"^ 

LcO(1->0)/'-FIR 



Fig.2. L[cii]/Lfir vs. Lco(i-0)/iFiR for LESS J033229.4 (red 
point with error bars dominated by the uncertainty in Lfir) com- 
pared with Galactic star-forming regions (crosses), starburst nu- 
clei (filled squares), non-starburst nuclei (open squares), normal 
galaxies (triangles), local ULIRGs (circles), and high-redshift 
sources (asterisks and diamonds). Objects known to have low 
metallicity are overplotted by red crosses. Black lines represent 
the solar metallicity PDR model calculations for gas de nsity (n) 
and FUV field strength (Go) from'Kaufman et alj jI999h . A vec- 
tor indicates how galaxy data points are shifted for the PDR anal- 
ysis in § 3. This figure is adapted from iStacev et al. (2010 |, with 
additio nal data on B RI 1335-0417 (z=4.4I) from .Riechers et all 
(l2008h and lWagg"et al. (2010), and on the low metallicity galaxy 
NGC 4142 (ICormier et al.ti201Q) . LESS J033229 is located near 
the high L[cii]/^co(i-0) region occupied only by low metalHcity 
galaxies. 



mal galaxies, and both nearby and l<z<2 starforming galaxies. 
However, it is significantly less intense than observed in Galactic 
OB star forming regions, three z>4 AGNs, and local ULIRGs 
(though the latter may also be influenced by dust-bo unded HII 
regions; lAbel et all 120091; iGracia-Carpio et al.L I20I1D . The un- 
certainties on the [CII]/FIR and CO/FIR ratios are large, but 
they are dominated by the uncertainty in Lfir, which may be 
1.5 times higher (Coppin et al., 2009). Such an increase in LpiR 
would shift the LESS J033229.4 point to the lower left, along 
the line of constant L[cii]/Lco(i-0)=10''. Similarly, an (unlikely) 
higher AGN contribution to Lfir would shift the point to the 
top right. Neither shift would change n or Go to a level that 
would change our conclusions. Our results are thus broadly 
consistent with the findings of IStacev et al.l (120101) who con- 
cluded that AGN dominated systems have ~8 times higher 
FUV radiation fields than starburst dominated systems. The 
n X! 10''cm"^ an d Go ~ 10-^ im ply a PDR surface temper- 
ature of -300 K (iKaufman et all Il999). Using formula 1 of 
IHailev-Dunsheath et alj ( 1201 Ol) . we estimate the atomic mass as- 
sociated with the C^ region Ma=1.0 +0.3x 10" ' Mq, i .e. similar 
to the M(H2)=1.6+0.3x 10'"Mo (Cop pin et al.[l2010l) . 

The very high L[cii]/Lco(i-0)~10'* places LESS J033229.4 
near the upper envelope of the solar metallicity plane-parallel 
slab PDR models in Figure |2l Only spherical PDR models 
can explain L[cii]/Lco(i-0)>10'' ratios, in particular wh en the 
metallicity in reduced (e.g. Bolatto et al., 1999; Rolli g et all 
|2006). Observational evidence is also building up that the high 
[CII]/CO region contains mainly low metallicity systems, where 
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the [CII] is enhanced relative to CO and the FIR dust contin- 
uum. Extreme L[cii]/ico(i-0) ratios up to 7x10"* have only been 
observed in low metallicity dwarf galaxies (e.g. Madden, 2000), 
and recently in star-forming region s of the low metallicit y irreg- 
ular galaxy NGC 4214 (see Fig. l2l ICormier et all I20T0I) . These 
high [CII]/CO ratios have been explained by the reduced dust 
abundance in low metallicity environments, which increases the 
mean free path of the FUV photons. This will not only decrease 
the FIR intensity and size of the CO cores, but also increase the 
size of the [CII] emitting regions as the photodissociating pho- 
tons traverse a larger volume of the molecular cloud. This ex- 
tended [CII] region may also contain self-shielded H2, which can 
boost the CO-to-H2 conversion factor by up to 1 00 times as see n 
in ICIO which has -1/5 solar metallicity (Ma dden et al.Lll997l) . 
The total gas mass based only on CO observations may thus be 
substantially underestimated in these low metallicity systems. In 
contrast, the atomic gas mass derived from the [CII] line only 
scales linearly with the abundance, and is less affected. While 
the L[cii]/Lco(\-u) ratio in LESS J033229.4 is not as extreme as 
in dwarf galaxies, such high ratios are still only observed in low 
metallicity environments. We therefore propose low metallicity 
as a natural explanation also in LESS J033229.4. It is also re- 
markable that three of the six z>2 [CII] detections have sim- 
ilarly high L,cni/Lcon-0) (F ig.|2ll, including BRI 0952-011^ 
where iMaiolino e t al.l (12009^ ^ already proposed low metallicity 
as an explanation for the large L[cii]/iFiR ratio. Unfortunately, 
we do not have independent measurements of the global metal- 
Ucity in LESS J033229.4 or any of the high-z objects with [CII] 
detections, though lower gas metallicities have indeed be en ob- 
sei-ved in z>3 star forming galaxies (Maiol ino et al.LlToOSb . If, as 
expected, the highest redshift galaxies indeed have lower metal- 
licities, their longer FUV mean free path lengths imply that for a 
constant FUV flux and density, the Go will be lower d ue to an in- 
creased geometrical dilution of the FIR flux (llsrael & Malonev, 
I2011h . The [CII]/FIR ratio will thus also increase (Fig. |2). We 
conclude that [CII] rather than CO should be the prime tracer of 
the ISM in the highest redshift (non AGN-dominated) galaxies, 
because [CII] will not only be brighter, but also provides a more 
reliable determination of the total gas mass. 

The high redshift star-formation dominated objects have 
similar Go fields than local star forming regions (Fig. |2), but 
much higher FIR luminosities. So the foimer appear to be scaled 
up versions of local starbursts. As the FUV emission is ab- 
sorbed by the dust in the molecular clouds, and re-radiated in 
the FIR, the extent of t he emitting r egion for a given Lfir will 
Stacev etal] (l2( 



increase for smaller Gi 



■2010.) use two scaling laws 

from lWolfire et al.l (1 19901) to constrain this size: Gq oc ALm/D^ 
if the mean free path of a FUV photon A is very small and 
Go Lm/D^ if the mean free path is very large. Applied to 
LESS J033229.4, and assuming Go~10^, this yields a physi- 
cal size DpDR=L6-3.7kpc. As argued above, in low metallic- 
ity environments, the larger FUV mean free path approximation 
would be more appropriate, so the actual value is probably at 
the higher end of this range. This scale size is remarkably sim- 
ilar to the 4kpc full width half maxi mum of the s t ellar popu- 
lation derived from /iT-band imaging dCoppin et all 12010 ). and 
shows that the starburst in LESS J033229.4 is affecting the en- 
tire galaxy. Such kpc-scale starburst have been found by several 
authors (e.g. lBiggs & Ivisonil2008HCarilli e t al.l. 1201 0|K_ with in- 
dividual components as small as lOOpc Hwinbank etalll2010l) . 



^ The lensing amplification in BRI 0952-0115 also implies a less 
luminous AGN, which may explain why BRI 0952-0115 is closer to 
the starburst rather than AGN dominated systems in Fig.|2] 



and sometimes pa r t of a merging system separated by a few 
kpc (iTacconi et al.L l2008h . However, the presence of a compact 
AGN, and differential lensing amplification of these different 
regions has complicated the interpretation. High spatial resolu- 
tion [CII] and FIR observations of unlensed starburst dominated 
systems with the Atacama Large Millimeter and submillimeter 
AiTay will provide a major progress in this field. 



4. Conclusions 

We have detected bright [CII] emission in a z=4.76 SMG, where 
the AGN contributes Uttle the FIR and [CII] luminosities. We 
find an atomic mass derived from [CII] comparable to the molec- 
ular mass derived from CO. By similarity to local dwarf galax- 
ies, the high L[cii]/Lco(i-0)~10'' ratio suggests a low metallicity 
in the highest redshift star-forming galaxies. The enhanced [CII] 
suggests that, as long as the AGN contribution is small, the [CII] 
line may be the most powerful and reliable tracer of the ISM. 
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